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Based on the analogy with superconductor physics we consider a scalar-vector-tensor gravita¬ 
tional model, in which the dark energy action is described by a gauge invariant electromagnetic 
type functional. By assuming that the ground state of the dark energy is in a form of a condensate 
with the U(l) symmetry spontaneously broken, the gauge invariant electromagnetic dark energy 
can be described in terms of the combination of a vector and of a scalar field (corresponding to 
the Goldstone boson), respectively. The gravitational field equations are obtained by also assuming 
the possibility of a non-minimal coupling between the cosmological mass current and the supercon¬ 
ducting dark energy. The cosmological implications of the dark energy model are investigated for a 
Friedmann-Robertson-Walker homogeneous and isotropic geometry for two particular choices of the 
electromagnetic type potential, corresponding to a pure electric type held, and to a pure magnetic 
held, respectively. The time evolutions of the scale factor, matter energy density and deceleration 
parameter are obtained for both cases, and it is shown that in the presence of the superconducting 
dark energy the Universe ends its evolution in an exponentially accelerating vacuum de Sitter state. 

By using the formalism of the irreversible thermodynamic processes for open systems we inter¬ 
pret the generalized conservation equations in the superconducting dark energy model as describing 
matter creation. The particle production rates, the creation pressure and the entropy evolution are 
explicitly obtained. 

PACS numbers: 03.75.Kk, 11.27.-|-d, 98.80.Cq, 04.20.-q, 04.25.D-, 95.35.-|-d 


I. INTRODUCTION 

The A Cold Dark Matter (ACDM) model of cosmology 
is remarkably successful in accounting for most of the ob¬ 
served features of the Universe. The recent Planck satel¬ 
lite data from the 2.7 degree Cosmic Microwave Back¬ 
ground full sky survey [l|, have generally confirmed 
again the present day standard cosmological paradigm. 
However, a number of fundamental questions at the very 
foundations of cosmology and gravitation still remain 
open, and unanswered. Perhaps the most important chal¬ 
lenge facing modern cosmology is the understanding of 
the mechanism of the acceleration of the late universe, 
which is usually attributed to the presence of the mys¬ 
terious dark energy. Q. In fact, as fundamental candi¬ 
dates responsible for the cosmic expansion, the standard 
ACDM model of cosmology has favored dark energy mod¬ 
els involving time-dependent scalar fields. Scalar fields 
naturally arise in many particle physics models, including 
string theory. On the other hand, the underlying dynam¬ 
ics of inflationary models, assumed to be of fundamental 
importance for the understanding of the early history of 
the Universe, also depend essentially on a single scalar 
field, the inflaton, rolling in some underlying potential 
0. The possibility that a single canonical scalar field cj), 
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with a non-zero potential, called quintessence, could be 
responsible for the late-time cosmic acceleration, was also 
explored in much detail The well-known action for 

a scalar held in the presence of gravity is 


S^ = 






( 1 ) 


where R is the Ricci scalar, G is the gravitational con¬ 
stant, and V{(f>) is the self-interaction potential, respec¬ 
tively 0. 

In opposition to the behavior of the cosmological con¬ 
stant, the quintessence equation of state changes dynam¬ 
ically with time In fact, many other exotic huids 
have been proposed to explain the accelerated expansion 
of the Universe. Some of the proposed models are the 
so-called fc—essence models, in which the late-time accel¬ 
eration is driven by the kinetic energy term of the scalar 
held 0- A number of coupled models, where dark en¬ 
ergy interacts both quantitatively and qualitatively with 
dark matter, have also been proposed [Q, as well as uni- 
hed models of dark matter and dark energy [l^ . For a 
review of the dark energy candidates see [1^. 

An intriguing alternative about the nature of dark en¬ 
ergy, which was also intensively investigated in the lit¬ 
erature, is the possibility that it could be described by 
a vector held, which can be at the origin of the present 
stage of cosmic acceleration. In its simplest formulation 
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the action for the vector field dark energy model is 


Sv= 


{ R 


16ttG 


-E 


+ V{A^) 


Lr. 


( 2 ) 


where = g^^’'AlAl, and L„ is 

the matter Lagrangian [l4| . This vector (or more exactly 
Yang-Mills) type action for the dark energy thus contains 
three identical components obtained by generalizing the 
Lagrangian of a single vector field. The term V{A^) is 
a self-interaction potential that explicitly violates gauge 
invariance. The cosmological implications of the vector 
type dark energy models have been investigated in [l5|. 

More general vector field dark energy models, in which 
the vector field is non-minimally coup led to the gravita¬ 
tional field, have been proposed in [1^. By assuming that 
the Universe is filled with a massive cosmological vector 
field, with mass /ta, which is characterized by a four- 
potential A^{x''), = 0,1,2,3, which couples non- 

minimally to gravity, and by introducing, in analogy with 
electrodynamics, the field tensor — Vi,A^, 

the action for the non-minimally coupled vector dark en¬ 
ergy theory can be written as 


5 = - 


R + + -AiiA^A'^ + wA^A'^i? + 


riA>^A''Rf,^ + 167rGoL„ 


y/^dn, 


( 3 ) 


where R^^^ is the Ricci tensor and Gq is the gravitational 
constant. In Eq. ([3]) uj and rj are dimensionless coupling 
parameters. 

At first sight the gravitational actions given by Eqs. (HD 
and (ED look totally different, from both mathematical 
point of view, as well as from the physical interpretation 
point of view. However, they can be in fact interpreted 
and understood as the limiting cases of a single physical 
model, related to the spontaneous breaking of the elec¬ 
tromagnetic U(l) symmetry. Thus an approach is used 
to describe superconductivity from a fundamental point 
of view (Trl - fl^ . 

Erom a general physical point of view in the theory 
of superconductivity the existence of a quantum con¬ 
densate (superconducting state) is described by a non- 
vanishin g va lue of a gauge dependent complex order pa¬ 
rameter |17Hl9l |. In bosonic systems superfluid behavior 
occurs when the expectation value of the bosonic field 
parameter tp has a nonzero value, (■0) ^ 0. On the other 
hand the existence of superconductivity is also induced 
by a nonzero value of the expectation value of the pair 
field operator. Therefore, in the ground state of a su¬ 
perconducting system a quantum condensate (eQ,/ 30 “ 0 ^) 
forms (TtI - Di^ . Since the difermion operator has charge 
—2e, the important result that the quantum condensate 


breaks the electromagnetic U(I) symmetry is obtained. 
Another fundamental quantity in the model is a scalar 
field, $, which plays the role of the order parameter. 
Under a gauge transformation A^^ ^ A^ + d^A, the 
scalar field transforms like the condensate wave func¬ 
tion 0 —>• => —>■ Note that in the 

zero temperature superconductivity theory one also in¬ 
troduces the Goldstone field 0 as the phase of the field 
d), d) = as well as the gauge invariant Fermi fields, 

0 = [13,El- 

Hence from a fundamental point of view a supercon¬ 
ducting system can be described by a gauge invariant 
Lagrangian, depending on the wave function 0, and on 
the vector potentials A^ and A simplified model is 

obtained after integrating out the Fermi fields. Thus one 
obtains a gauge invariant Lagrangian, depending only 
on and respectively. The important require¬ 

ment of the gauge invariance of the theory implies that 
these bosonic fields must appear only in the combinations 
= A/^A^ — A/vAf^ and —V^0, respectively. There¬ 
fore the Lagrangian describing a superconductor from a 
fundamental physical point of view has the form [13,13 

L = -\J F^,F^^^d^r + Ls (A^ - V^0), (4) 

where Lg (A^ — V^0) is an arbitrary function of the argu¬ 
ment A^ — The only physical condition required on 
the superconductor Lagrangian Ls is that in the absence 
of A'^ and 0 it gives rise to a stable state of the system. 
In particular, this requires that the point A^ = is 
a local minimum of the theory (this property can fully 
explain the Meissner effect in superconductivity theory 
[l7|b Therefore, we require that the second derivative of 
the superconductor Lagrangian Lg with respect to its ar¬ 
gument must be nonzero at the point A^ = Em. 

It is the goal of the present paper to consider a gravi¬ 
tational model in which dark energy is described by a La¬ 
grangian of the form given by Eq. ED, resulting from the 
breaking of the U(I) symmetry in the ground state dark 
energy condensate. By analogy with condensed mat¬ 
ter physics we call this model the superconducting dark 
energy model. The gravitational field equations of the 
model are derived from an action principle, and the cos¬ 
mological implications are investigated in a background 
homogeneous and isotropic flat Friedmann-Robertson- 
Walker geometry. We consider two distinct classes of cos¬ 
mological models, corresponding to two different choices 
of the electromagnetic potential A^ of the dark energy. 
In the first model A^ has only a non-vanishing tempo¬ 
ral component, while in the second case we assume non¬ 
vanishing spatial (magnetic) components of the poten¬ 
tial. In both cases we assume that the dark energy self¬ 
interaction potential is constant. 

In the present dark energy model, due to the coupling 
between the matter current and the electromagnetic and 
scalar potentials of the dark energy, the matter energy- 
momentum tensor is not conserved. By using the formal¬ 
ism of the open thermodynamic systems introduced in 
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(see also [1^ for recent investigations of particle 
creation in cosmology), we interpret the generalized con¬ 
servation equations in the superconducting dark energy 
model from a thermodynamic point of view as describ¬ 
ing irreversible matter creation processes. Thus in the 
present model particle creation corresponds to an irre¬ 
versible energy flow from the superconducting dark en¬ 
ergy to the created matter constituents (both normal and 
dark). We explicitly obtain the equivalent particle num¬ 
ber creation rates, the creation pressure and the entropy 
production rates. The temperature evolution laws of the 
newly created particles are explicitly derived. We also 
show that due to the superconducting dark energy - mat¬ 
ter current coupling, during the cosmological evolution a 
large amount of comoving entropy could be produced. 

The present paper is organized as follows. In Sec¬ 
tion El the gravitational field equations of the supercon¬ 
ducting dark energy model, a scalar-vector-tensor theory 
with broken U(l) symmetry, are derived from a varia¬ 
tional principle. The equations of motion of the scalar 
and vector fields are also obtained. The cosmological 
applications of the theory are investigated in Section Hill 
Two distinct dark energy models are considered: an elec¬ 
tric type, in which the vector potential has only a time 
component, and a magnetic type, with the vector po¬ 
tential having only spatial components. The cosmologi¬ 
cal properties of both models are investigated in detail. 
The thermodynamic interpretation of the superconduct¬ 
ing dark energy model is considered, in the framework of 
the thermodynamic of open systems and irreversible pro¬ 
cesses, in SectionllVl We discuss and conclude our results 
in Section |Vl In this paper we adopt the Landau-Lifshitz 
metric conventions, and we use the natural system 
of units with SttG = c = 1. 


II. FIELD EQUATIONS OF THE 
SUPERCONDUCTING DARK ENERGY MODEL 


In the following we assume that the interaction of the 
gravitational and of the superconducting dark energy 
scalar-vector fields is described by a Lagrangian which is 
required to satisfy the following standard conditions: a) 
the Lagrangian density is a four-scalar b) the free-held 
energies are positive-dehnite for all the metric, scalar 
and vector helds c) the resulting theory is metric and 
d) the held equations contain no higher than second or¬ 
der derivatives of the helds [2l|. Based on the analogy 
with superconductor physics we consider a gravitational 
scalar-vector-tensor action of the form 


S = - 


R 

'2 




{A^ - V^</.) {A, - V,<).) + U (A^, (/.) - 


Lm {gp.u, ■*/') 




where A and a are constants, Lm{gfiiy,tp) is the La¬ 
grangian of the total (ordinary baryonic plus dark) mat¬ 
ter, and = pu^ is the total mass current, where p is 
the total matter density (including dark matter), and 
is the matter four-velocity. We assume that the bary¬ 
onic and dark matter are comoving. The third term in 
the action Eq. ([5]) follows from the assumption that the 
superconducting dark energy is close to the minimum 
A/i = In this case the general superconducting La¬ 

grangian O can be expanded in power series as [13,13 


Ls {Af^ — V^(/)) ~ Tq+x 




2 5{A^-y^<l>y 


{Afj, — yficj)) -I-..., 


( 6 ) 

where Lq is a constant. Hence the superconducting type 
Lagrangian Lg (A^ — V^^) gives a quadratic contribu¬ 
tion in Afj_ — V to the gravitational Lagrangian. We 
have also assumed the possibility of an interaction be¬ 
tween the total matter flux and the superconduct¬ 
ing dark energy gauge invariant potentials A^ — Vfji,4>. 
V{A^,(f)) is the self-interaction potential of the scalar 
and vector fields, with A^ = A^A^, in which we have 
also included the constant Lq. When 0 = 0, that is, 
the scalar field vanishes, the action ([5]) gives the pure 
vector model of the dark energy. When the electromag¬ 
netic type potential A^ = 0, we recover the standard ac¬ 
tion of the minimally coupled scalar-tensor theory. Hence 
the gravitational action ([5]) gives a unified framework for 
the minimal inclusion into the gravitational action of the 
scalar-vector interactions, under the assumption of the 
existence of a U(l) broken symmetry. The second and 
third terms in the gravitational action are also similar to 
the Stueckelberg Lagrangian [T7I| . 

We define the energy-momentum tensor of the matter 


as 


T 


2 ' d{y/^Lm) 


d d{yf^Lm) 
dx^ d /dx^) 


■ (7) 


By making the important assumption that the La¬ 
grangian density Lm of the matter depends only on the 
metric tensor components and not on its derivatives, 
we obtain the expression = Lmg^v — ‘^dLm/dg^''. 

By varying the action ([5]) with respect to the metric 
tensor we obtain the gravitational field equations for the 
superconducting dark energy model as 


Rfii^ — ■^Rg^^u — 701/ + ^ ^ —-I- -FapF°‘dg^^'^ _|_ 

A (A^ - V^0) (A, - Vi/0) - ^ (Al“ - V“0) (A„ - V„0) X 

ffAii/ + ajfj, {A^ - Vi.0) - {Ap - S/pcj)) -h 

V{A^y>)g^,, (8) 


where T^j/, the energy-momentum tensor of the ordinary 
matter, is given by 


701 / = {p + p) - pgfj.^, 


(9) 
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where p is the total thermodynamic pressure of the mat¬ 
ter components (baryonic and dark). By taking the vari¬ 
ation of the action Eq. (I5|) with respect to the scalar field 
(j) we obtain 


Sd,S = - 




a . 




^—g(Kl. (10) 


With the use of the mathematical identity 

it follows that if the conditions 


( 11 ) 


V^^'^ = 0,V^j'^ = 0, (12) 

are imposed on the dark energy vector potential, and 
on the baryonic matter flow, the variation in the integral 
m of the first term, containing A^, and of the last term, 
containing j^, vanish identically. Therefore we obtain 
the result that if the conditions given by Eqs. m are 
satisfied, then the scalar field satisfies the standard Klein- 
Gordon equation, 

+ d^V {A^ (f) = 0. (13) 

This case corresponds to the minimal coupling of the 
scalar and vector fields. However, in the following we will 
use a more general approach, in which no additional con¬ 
straints are imposed on the fields or on the hydrodynamic 
flow. Therefore, the variation of the action Eq. gives 
the following coupled evolution equation for the scalar 
and vector fields, 

Xg^^''V^V,c^+d^V{A\4>)-XV^A^^-^V^j^^ = 0. (14) 

By varying the superconducting dark energy action 
Eq. ® with respect to we obtain first 


5ai^S = 


_ - Xg>^’' {A^ - y^cj}) SA^ - 

47r 


-j>^SA^ + 2dA2V{A^,cf) Af^SA^ 


y/—gdVL = 0. 

(15) 


By taking into account the identity Vy (F^'^dH^) = 
X/yF^'^SA^ + ySAfj_, after partial integration and 

the use of Gauss’ theorem, it follows that the supercon¬ 
ducting dark energy vector field satisfies the equation 

(16) 

47r 

where 


J'" = 


{Ay 


V.</)) + 


2dA^V {A^,4>)A>^ . 

(17) 


The divergence of the dark energy field tensor can be 
obtained as yyF^^’^ = {l/^/^) dy {y/-^F^''). 

By its definition the dark energy electromagnetic type 
tensor F^^'^ satisfies the Bianchi identity 

pF^y = 0, (18) 

where ^3 ^;]^e complete antisymmetric unit tensor 

of rank four. 

Finally, by taking the covariant derivative of the field 
equations Eqs. ® we obtain the matter conservation 
equation in the presence of a superconducting dark en¬ 
ergy as 

ir (Ay - yy<p)] - ^yy/ (Af, - + 

d^V (H2, (f) + 2dA^V (H2, A^y^Ay = 0. (19) 

The derivation of Eq. (ITO)) is presented in Appendix 
By taking into account the explicit form of the energy- 
momentum tensor, given by Eq. we obtain 

(-I- y^p) u^Uy + {p + p) Uyy^u^ + {p + p) u^y^Uy - 

(A, - yy^)] - ^yyf (Ap - Vp0) + 

d^V (A2, <(,) A, + 2dA2V A“V„A, = 0. (20) 

By multiplying Eq. dsni) with and by taking into ac¬ 
count the mathematical identity u’^y^Uy = 0 we obtain 
the energy conservation equation in the superconducting 
dark energy model as 

p + 3{p + p)H + ^u’^y^ \j>^ (Ay - V,0)] - 

f S [/ (^/3 - ^0^)] + u^'Ay + 

2dA^V {A^,(j))u''A^yaAy=Q, ( 21 ) 

where we have introduced the Hubble function iJ = 
(l/3)V'^u^, and we have denoted ' = = d/ds, re¬ 

spectively, where ds is the line element corresponding to 
the metric g^y, ds^ = g^ydx^dx'' . 

By multiplying Eq. (EU)) with the projection operator 
defined as /i)( = S'^ — uxu'', and satisfying the relation 
Uyh'^ = 0 , gives the momentum balance equation for a 
perfect fluid in the superconducting dark energy model 
as 

A li/A f 

ir (Ay - yy^)] + ^yyf (A^ - V^^) - 

d^V (A 2 , <^) A, - 2dA^V (A 2 , A“V„A,|. ( 22 ) 

III. COSMOLOGICAL APPLICATIONS 

We assume that the metric of the Universe is given 
by the isotropic and homogeneous Friedmann-Robertson- 
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Walker metric, 

ds^ = dt^ — a^{t) {dx^ + dy^ + dz^) , 


As an independent variable we introduce, instead of the 
cosmological time t, the redshift z, defined as 1 + z = 1/a. 
(23) Therefore 


where a{t) is the scale factor describing the expansion of 
the Universe. We assume that the cosmological matter 
is comoving with the cosmological expansion, and there¬ 
fore we choose the four velocity of the cosmological fluid 
as = (1,0,0, 0). Hence the components of the four- 
current vector are = (p, 0,0,0). In the Friedmann- 
Robertson-Walker geometry the Hubble function takes 
the form H = a/a, since = ' = d/dt. To describe 

the decelerating/accelerating nature of the cosmological 
expansion, we use the deceleration parameter q, with the 
definition 


dH dH dz ^ 

dt dz dt dz 


(30) 


As a function of the redshift the deceleration parameter 
is obtained as 


q = {l + z) 


1 dH{z) 
H{z) dz 


- 1 . 


(31) 


In the following we will explicitly investigate two dis¬ 
tinct superconducting dark energy models. 
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d I H 

dlH~ ~~1P 


(24) 


A. Electric dark energy models 


Moreover, from the homogeneity of the Universe it 
follows that the scalar and vector fields </> and can 
be only functions of the cosmological time t, so that 
(j) = (j){t) and Afj, = A^{t) = {Ao{t),Ai{t),A 2 {t),A 3 {t)), 
respectively. The non-zero components of the dark en¬ 
ergy tensor are given by F^it) = -Ai{t), i=l,2, 3, 
and F*°(t) = Ai{t)/a?{t), i = 1,2,3. Hence we obtain 


FayF°^^ = — (2/a^(t)) ^i(^) ■ Then, the cosmo¬ 

logical equations corresponding to the superconducting 
dark energy model are given by 


3^2 



1 

87ra^(t) 



A 

2 



A 

2a? it) 



+ V{A\<t>). 


+ 

(25) 


We assume that the dark energy vector potential has 
the form = (Ao(t), 0, 0,0), that is, the dark energy 
vector potential has only one, electric type, component. 
For this choice = 0, V/r, v G [0,1, 2, 3]. The gravita¬ 
tional field equations describing the cosmological dynam¬ 
ics in the presence of the superconducting dark energy 
take the form 





+y {A‘^,(j)j , 

(32) 


2H+3H^ = -p- 




(33) 


^ |a^ A (Ao - 0 ) + I - a^d^V (A^ </) = 0, (34) 


2H + = -p 

A? 


A 


aHt) 

a 


(Ao- 




1 1 
47ra^(t) 87ra^(t) 




\k^l 


X 


2a‘^{t) 


\k^l 


(Ao-</)p + U(A2,0),f = 1,2,3, 


(26) 


A (Ao - </) + |p - 2dA2V (A2,0) Ao = 0. (35) 

In order to close the system of equations (l32])-([35|) the 
baryonic equation of state p = p (p) must also be pro¬ 
vided. In the following we will restrict our analysis to the 
case of a constant self-interaction potential of the super¬ 
conducting dark energy field, V (A^, </) = Vq = constant. 
Then from Eqs. (l34l) and (I35|) we obtain 


A(/- 3 


p + A 


^Aq — (fj 


H-\A3--p- 


a^U(A2,</)) =0, 


A ^Aq — (/)^ + —p — 29^2U (A^, (/) Aq — 0, 


(27) 

(28) 


,Ao-^]=--p. 


(36) 


Hence the generalized Friedmann equations of the cos¬ 
mological expansion in the presence of the electric type 
superconducting dark energy become 


— p — + To — p + Pde, 


(37) 


Afc + FAfe(t) + 47rAAfc(t) — 87r9yi2U (A^, (/) Afc(t) — 0, _ ^2 

fc = l,2,3. (29) 2H + 3H^ = -p+ —p^ + Vo = -P-PDE, (38) 
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where we have denoted 


PDE = -^P 


Vo, 


By introducing a set of dimensionless variables 
{9,T,h,vo), defined as p = (8A/a^) 0, t = a/V24Ar, 

(39) H = h, Vo = (a^/8A) Vo, Eq. (H51) becomes 


and 

2 

Pde = ~'^P^ ~ ^0, (40) 

respectively. From Eqs. dSzl) and (l38l) we obtain 


2H = -{p + p) + -p^. 


(41) 


The energy conservation equation can be written as 


dt 




(42) 


The deceleration parameter can be obtained as 

_ (P + 3p) - (a!^/2A) p^ - 2Vo 
2[p-{ay8X)p^ + Vo] ■ 


For a dust Universe with p = 0, the deceleration pa¬ 
rameter takes the form 


p - (q;^/2A) p^ - 2Vo 
2 p - (aV8A) 


From the field equations Eqs. (l^ - (l4T|l we obtain the 
time evolution equation of the baryonic matter density 
as 


._ v'P - (aV8A) p"^ + Vo{p + p[l- (aV4A) p] ] 
[l-(aV4A)p] 

(45) 

In terms of the redshift z, the density evolution equa¬ 
tion of the baryon density in the electric type supercon¬ 
ducting dark energy models is given by 

dp{z) ^ 3 p{z) + p{z) [1 - (aV4A) p{z)\ 

dz 1 + z 1 — (a^/4A) p(z) 


In order to characterize the dark energy, and its evolu¬ 
tion properties, we also introduce the dark energy equa¬ 
tion of state parameter wde, defined as 


wde 


Pde 

Pde 


- (a^/8A) -I- Vo 

(aV8A) p-^ + Vo 


(47) 


1. Dust Universes with electric type superconducting dark 
energy 

In the case of the dust Universe, with p = 0, Eq. (1451) . 
describing the time dynamics of the matter density in 
the presence of the superconducting electric type dark 
energy takes the form 

p = -'/Spy^p- (Q;2/8A)p2 + Vo. (48) 


^ = -eVd-e^+vo, (49) 

dr 

while the dimensionless Hubble parameter h is obtained 
as 


h = yl0-e^+vo. (50) 

In the dimensionless time variable r we have /i(t) = 
[3/a(T)] {da/dr). The deceleration parameter of this 
model is given by 


1 0 - 46»2 - 2vo 

2 e-0‘^ + vo 


(51) 


while the dark energy equation of state parameter wde 
can be obtained as 


Wde = — 


- 0 ^ + ^0 
6^ + Vo 


(52) 


The time variation of the redshift z can be obtained 
from the equation 


dz .Ida 1 + z r- - — - 

— = -{\ + z)-— = - —\J0- 0^ +vo. (53 

dr a dr 3 

In terms of the redshift z the evolution of the dust 
electric type superconducting dark energy Universe is de¬ 
scribed by the simple relations 

p{z) = po{l + zf , (54) 


Vs [ 


2 2 
« Po, 


Po(l + 2)^-^(l + ^)^ + Uo 


1/2 


(55) 


^(^)^l„o(l+.)^-(aVg/2A)(l + ./-2r„_ (56) 

2 po(l + zf - {o^pI/SX) (1 -b zf + Vo 

where po is the matter density of the Universe at the 
present time z = 0. 

The variations with respect to the redshift z of the 
Hubble function h of the Universe, of the matter energy 
density 0, of the deceleration parameter q, and of the 
parameter of the dark energy equation of state are rep¬ 
resented, for different values of vo, in Figs. [TJ41 The 
initial values used to numerically integrate the cosmolog¬ 
ical evolution equations are 0(0) = 0.1, z(0) = 5, and 
a(0) = 1/6, respectively. 

As one can see from Fig. [1] the Hubble function h of 
the Universe is a monotonically increasing function of the 
redshift (monotonically time decreasing function). In the 
early stages of evolution, at around z « 5, h is basically 
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FIG. 1: Variation with respect to the redshift z of the dimen¬ 
sionless Hubble function h of the electric type superconduct¬ 
ing dark energy filled Universe with p — 0 for different values 
of Vo’- Vo = 0.001 (solid curve), vo = 0.002 (dotted curve), 
Vo = 0.003 (short dashed curve), vo = 0.004 (dashed curve), 
and Vo = 0.005 (long dashed curve), respectively. 
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FIG. 2: Variation with respect to the redshift z of the dimen¬ 
sionless matter energy density 0 of the electric type supercon¬ 
ducting dark energy filled Universe with p = 0 for different 
values of no: vo = 0.001 (solid curve), vo = 0.002 (dotted 
curve), Vo = 0.003 (short dashed curve), vo — 0.004 (dashed 
curve), and vo = 0.005 (long dashed curve), respectively. 
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FIG. 3: Variation with respect to the redshift z of the de¬ 
celeration parameter q of the electric type superconducting 
dark energy filled Universe with p = 0 for different values 
of Vo’. Vo = 0.001 (solid curve), vo — 0.002 (dotted curve), 
Vo = 0.003 (short dashed curve), vo = 0.004 (dashed curve), 
and Vo = 0.005 (long dashed curve), respectively. 
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FIG. 4: Variation with respect to the redshift 2 of the param¬ 
eter wde of the dark energy equation of state of the electric 
type superconducting dark energy filled Universe with p — 0 
for different values of vo’. uo = 0.001 (solid curve), vo = 0.002 
(dotted curve), vo = 0.003 (short dashed curve), vo = 0.004 
(dashed curve), and vo = 0.005 (long dashed curve), respec¬ 
tively. 


independent on the values of vq. The matter density of 
the Universe, presented in Fig. [2 is a monotonically in¬ 
creasing function of z, tending in the small z limit to zero, 
lim 2 _j.o 0(z) = 0. Its evolution is basically independent of 
the range if the numerical values of vq ■ The redshift vari¬ 
ation of the deceleration parameter q, depicted in Fig. [3 
shows that in the present model the Universe starts from 
a decelerating phase at around 2 ~ 5, with q having 
values of the order of g ~ 0.2 — 0.3. This initial value 
increases in the early stages of the cosmological evolu¬ 
tion, showing a decelerating expansion. At z ~ 1 — 2, 
the Universe starts to accelerate, with the decelerating 
parameter slightly decreasing and taking negative vales 
q < 0. The values of the deceleration parameter gradu¬ 
ally decrease with decreasing z, and the Universe enters 
into an accelerating stage, ending its evolution in a de 
Sitter stage, with gs;—latz~0. The parameter wde 


of the equation of state of the dark energy, presented in 
Fig. m starts with positive values, and, with decreasing 
z, it takes negative values. In the small redshift limit it 
tends to the value wde = — 1- 

The present day numerical values of the Hubble func¬ 
tion Hq and of the deceleration parameter go can be ob¬ 
tained as 


and 


Ho = 


[po - (aVo/8A) + Vb] 
V3 


1 [po - (q^P§/ 2A) - 2Uo] 

2 [po - (q:^P§/ 8A) + Uo] 


(57) 


(58) 


respectively. Therefore the free parameters a and A of 












the superconducting electric type dark energy model can 
be obtained from astronomical observations. 

Eq. (1491) can also be solved exactly, and thus we obtain 
the density as a function of time in an exact analytical 
form as given by 


0{t) = 


^\]Vo (- 6*0 + ^0 + Vo) + 


(00 + 2uo)e"'>v^ 


[9l + 8^2 - 4(00 - 2)0o^;o] x 


g2ro_ 45)gerV^Y^uo {-e^ + 00 + Vo) + 



el{Avo + + 4(00 + x 

\JVo (-^0 +^o+ Vo) - 

20o(0o + 2uo)ev^("+"«)| , (59) 


FIG. 5: Variation with respect to the redshift 2 of the Hub¬ 
ble function h of the radiation fluid Universe with p = pjS 
in the presence of electric type superconducting dark energy 
for different values of vg: vo = 0.001 (solid curve), uo = 0.002 
(dotted curve), vo = 0.003 (short dashed curve), vo = 0.004 
(dashed curve), and vg = 0.005 (long dashed curve), respec¬ 
tively. 


where we have used the initial condition 0 (tq) = 00 - 


The radiation fluid Universe with electric type 
supercondueting dark energy 


For a high density radiation fluid Universe, with matter 
equation of state satisfying the condition p = p/3, in 
the presence of electric type superconducting dark energy 
the basic evolution equation of the dimensionless matter 
density 0 is given by 


de 

dr 


t o . 2 , g/3 + 0(l-20) 

( 1 - 20 ) • 


(60) 


The deceleration parameter of the radiation Universe can 
be obtained as 



z 


20 - 402 - 2z;o 
2 (0 - 02 + no) ■ 


(61) 


The variations with respect to the redshift 2 of the 
Hubble function, of the energy density, of the decelera¬ 
tion parameter and of the parameter of the dark energy 
equation of state of the radiation fluid Universe in the 
presence of the electric type superconducting dark en¬ 
ergy are presented in Figs. [S][5] The initial conditions 
used to numerically integrate the cosmological evolution 
equation are 0(0) = 0.45 and 2 ( 0 ) = 25. 

We assume that the Universe was radiation dominated 
in the redshift range 5 < 2 < 25. The dimensionless Hub¬ 
ble function h of the high density Universe, presented in 
Fig. m is a monotonically increasing function of the red¬ 
shift (monotonically decreasing in time), while the energy 
density, depicted in Fig. [51 increases monotonically with 
the redshift 2 during the cosmological evolution. The de¬ 
celeration parameter, shown in Fig. [71 has positive values 
for the redshift interval 9 < 2 < 25, indicating a deceler¬ 
ating expansion. For large values of vq the deceleration 


FIG. 6: Variation with respect to the redshift 2 of the matter 
energy density 0 of the radiation fluid Universe with p — p/3 
in the presence of electric type superconducting dark energy 
for different values of vg: vg = 0.001 (solid curve), vg = 0.002 
(dotted curve), vg = 0.003 (short dashed curve), vg = 0.004 
(dashed curve), and vg = 0.005 (long dashed curve), respec¬ 
tively. 


parameter can reach the zero value at redshifts as high 
as 2 « 9, limz_).g 905 ~ 0. The time variation of 
the cosmological parameters h and 0 is practically in¬ 
dependent of the adopted small values of the parameter 
Vq. The parameter wde of the dark energy equation 
of state is represented in Fig. For large redshift val¬ 
ues 15 < 2 < 25, Wde is positive, while for 2 ~ 5 it 
approaches the value wde = ~ 1 , showing that in the 
present model the Universe becomes dark energy domi¬ 
nated at around 2 Ri 5. 
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FIG. 7: Variation with respect to the redshift 2 : of the deceler¬ 
ation parameter q of the radiation fluid Universe with p = pj?) 
in the presence of electric type superconducting dark energy 
for different values of uq: wo = 0.001 (solid curve), vo — 0.002 
(dotted curve), vg = 0.003 (short dashed curve), vo = 0.004 
(dashed curve), and vg = 0.005 (long dashed curve), respec¬ 
tively. 



FIG. 8: Variation with respect to the redshift 2 : of the pa¬ 
rameter of the dark energy equation of state wns of the radi¬ 
ation fluid Universe with p = p/3 in the presence of elec¬ 
tric type superconducting dark energy for different values 
of vg: vg — 0.001 (solid curve), vg = 0.002 (dotted curve), 
vg = 0.003 (short dashed curve), vg = 0.004 (dashed curve), 
and Vg = 0.005 (long dashed curve), respectively. 


S. The unified picture of the evolution of the Universe in 
the electric type superconducting dark energy model 

Finally, to conclude the investigation of the electric 
type superconducting dark energy model, we present a 
unified picture of the evolution of the Universe for the 
redshift range 2 G [0,25]. The variations of the Hubble 
function, matter energy density, deceleration parameter, 
and the parameter of the dark energy equation of state 
are plotted in Figs. 

To study the evolution of the electric type supercon¬ 
ducting dark energy we adopt for the redshift 2 the range 
from 0 to 25. We assume that in the range 2 € [5, 25] the 
matter content of the Universe can be (at least approx¬ 
imately) described by a radiation type equation of state 


FIG. 9: Variation with respect to the redshift 2 £ [0, 25] 
of the Hubble function h of the Universe filled by electric 
type superconducting dark energy, for different values of uq: 
Vg = 0.001 (solid curve), vg = 0.002 (dotted curve), vg = 
0.003 (short dashed curve), vg = 0.004 (dashed curve), and 
Vg — 0.005 (long dashed curve), respectively. 



FIG. 10: Variation with respect to the redshift 2 £ [0, 25] 
of the matter energy density 9 of the Universe filled with 
electric type superconducting dark energy, for different values 
of Vg-. Vg = 0.001 (solid curve), vg = 0.002 (dotted curve), 
Vg = 0.003 (short dashed curve), vg — 0.004 (dashed curve), 
and Vg = 0.005 (long dashed curve), respectively. 


p = pjZ. At 2 = 5 the Universe enters in the matter 
dominated era, with p 0. In this simplified model the 
transition from the radiation dominated era to the matter 
dominated phase is smooth, with all physical and ther¬ 
modynamical quantities continue at the transition point. 
Therefore the Hubble function and the matter density, 
represented in Figs. Inland 1101 are monotonically increas¬ 
ing functions of the redshift for the entire period. The 
deceleration parameter, shown in Fig. 1111 has a complex 
behavior. The Universe starts at 2 = 25 with a decel¬ 
eration parameter having a value of q ~ 0.10, and its 
early evolution is strongly decelerating, with q reaching 
the value g ~ 0.9 at 2 ~ 12, for vq = 0.001. For larger 
values of vq the expansion of the Universe is faster, with q 
reaching the value g = 0.6at2Ril4. After q has reached 
its maximum, it starts to decrease with decreasing 2 , and, 
depending on the numerical value of vq, reaches the value 
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FIG. 11: Variation with respect to the redshift 2 ; G [0, 25] 
of the deceleration parameter q of the Universe filled with 
electric type superconducting dark energy, for different values 
of Vo’. Vo = 0.001 (solid curve), vo = 0.002 (dotted curve), 
Vo = 0.003 (short dashed curve), vo = 0.004 (dashed curve), 
and Vo = 0.005 (long dashed curve), respectively. 



FIG. 12: Variation with respect to the redshift 2 G [0, 25] 
of the parameter of the dark energy equation of state wne 
for an Universe filled with electric type superconducting dark 
energy for different values of uq: vo = 0.001 (solid curve), vo = 
0.002 (dotted curve), vo = 0.003 (short dashed curve), vo = 
0.004 (dashed curve), and no = 0.005 (long dashed curve), 
respectively. 


q = 0 for 2 ~ 5 — 7. Then the deceleration parameter 
enters the negative range, with the Universe starting to 
accelerate at a higher rate, and reaching the value g Ri — 1 
(the de Sitter phase) at around 2 = 0. The parameter 
of the dark energy equation of state wbe, represented in 
Fig. 1121 is slowly decreasing from its maximum value 1 in 
the redshift range 2 ~ 17—25. For 2 < 15, wde decreases 
rapidly with decreasing 2 , and, depending on the numer¬ 
ical value of uo, reaches the limiting value wde ~ — 1 at 
2 « 7.5-10. 


4- Electric type superconducting dark energy Universe with 
conserved electric field and matter current 

Finally, we investigate the case in which the electric 
potential and the matter current satisfy the conserva¬ 
tion equations given by Eqs. (1121) . The time variation of 
can be immediately obtained from the Lorentz gauge 
equation imposed on which gives 

( 62 ) 

and 

A°{t) = (63) 

respectively, where Co is an arbitrary integration con¬ 
stant. The continuity equation of the matter hydrody¬ 
namic flow, {pu^) = 0 gives a similar dependence for 
the matter density p, 

P{t) = (64) 

where po is an arbitrary integration constant. The evo¬ 
lution of the scalar field (j) is decoupled from the electric 
and matter component, and in the presence of a constant 
potential, V = constant, follows a similar law as 

the electric and the matter fields, 

Ht) = (65) 

where fiio is an arbitrary integration constant. In the 
large time limit, all these fields tend to zero, and the 
Universe enters in an exponential, de Sitter type expan¬ 
sionary phase. However, the presence of the electric type 
superconducting dark energy modifies the cosmological 
dynamics of the Universe before it enters in the de Sitter 
stage. 


B. Magnetic dark energy models 


As a second superconducting vector type dark energy 
model we consider the case in which dark energy has a 
magnetic type structure, with its vector potential given 
by = ( 0 , Ai(f), A2(t), ^3(1)). In order to have an 
isotropic expansion the components of the superconduct¬ 
ing magnetic vector potential must satisfy the condition 
Ai{t) = A2{t) = ^3(1) = A{t). Hence for this choice 
of we obtain F)o = —Ai, = Aija}, i = 1,2,3, 
( 1 / 16 ^) = -( 3 / 8 ^)AVa 2 , ( 1 / 4 ^) ^0^0“ = 

— !^tt) A^!, and (l/47r)F^F“ = — (l/47r)A/a^ (no 

summation upon the index i). Therefore the gravita¬ 
tional field equations describing the isotropic and homo¬ 
geneous Universe in the presence of superconducting dark 
energy take the form 


3i7^ = 1 - 


3 A^ 
Stt a? 


( 66 ) 
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2 ct ; • 2 1 X 

2H-\-3H = —— ~4> — - -T " 1 “ ■ 


x^--p = o, 


+ K), (67) 


( 68 ) 


Therefore the system of gravitational field equations de¬ 
scribing the superconducting magnetic type cosmological 
dark energy model takes the form 


3 — — P - ^- 2 

8A 87r 


3 


+ K), 


(71) 


—- (aA) = —47rAA, 
a dt \ / 


(69) 


where for simplicity we have adopted a constant value Vq 
for the self-interaction potential of the scalar and vector 
fields, V = Vq = constant. The energy conserva¬ 

tion equation takes the form 


3 ( -b Vo ) a^d. 


(70) 


1 i 2 AA 2 


„d a? c? ; ^ 

2 -H—~o-(^ 2 ) 

a 8A 87r 2 a‘‘ 


A + -A + AirXA = 0. 

a 


From Eqs. (ED and ED we obtain 


2H = -{p + p) + -p^----X-. 


4A 


27r 


(73) 


(74) 


With the use of Eq. (1551) we can substitute the deriva¬ 
tive of the scalar field in terms of the matter density p. 


The deceleration parameter of the magnetic type super¬ 
conducting dark energy model can be represented as 


J 


9 = 


p + 3p— (a^/2A) + (3/47r) (^A^/a?-^ — 2Vo 


p - ((aV8A) p2) + (3/8^) + (3A/2) {A^/a?) + Eq 


(75) 


r 


while the parameter of the equation of state of the dark 
energy is given by 


wde = 


— (a^/8A) p^ + (3/87r) (^A ^+ (3A/2) jc?') -b Vo 

(a2/8A) _ (1/8^) (i2/a2) + (A/2) (Al2/a2) + Eq 

(76) 

In the case of a dust Universe, with p = 0, from the 
conservation equation Eq. EOD we obtain for the time 
derivative of the energy density p the equation 


P = 


^^ + |A^-^p 2 + 3p-3Eo 


Stt 3 ' 




H. 


(77) 


H = h,Vo = ( 8 A/a^) E)|, and by denoting 

a = Eo? 12, the system of Eqs. dTll - dTS]) can be writ¬ 
ten in a dimensionless form as 

1 /dE 


3h^ = 0 - -b 


dr 


Y? 

+ CT — 2 " + '^^Oi (78) 


dh 


2— —b 3h^ — 6^ — —— ( —) + —— 5 - + Vq, (79) 


^dh 
" dr 


= 0(20-1)--- 


1 1 fdYy crS^ 

3 \ dr y 3 

4 1 /dEA^ 2 E2 


3 a? \dT ) 3 a2 ’ 


(80) 


1. Dust Universes with magnetic type superconducting dark 

energy 


d^E , dE ^ „ 

, o—h h— -h crX j — 0, 

ar^ dr 


(81) 


For the case of dust matter, with negligible ther¬ 
modynamic pressure, we can take p = 0 in the 

gravitational field equations. Then by introduc¬ 
ing a set of dimensionless variables (d, E, r,/i, wq) de¬ 
fined as |p = ( 8 A/a^) 0,A= Y^ 87 r/ 3 E, t = r. 


dz , , 1 ^ 1 fdYY E2 

(82) 

The deceleration parameter and the parameter of the 
equation of state of the dark energy of the Universe filled 
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with magnetic type superconducting dark energy are ob¬ 
tained as 


0 _ 402 + 2E'Va' - 2no 

2(0-02+S'2/a2+^S2/a2+^^^)’ 


-02 + S'2/a2 -h 0-S2/a2 -h r;o 
02-S'2/3a2-ho-S2/3a2-huo’ 


(84) 


where a prime denotes the derivative with respect to the 
dimensionless time r. By taking the derivative with re¬ 
spect to T of Eq. dill), and with the use of Eqs. (1801) 
and m we obtain for the time variation of the matter 
density the equation 


where Oq is an arbitrary constant of integration. 

The variations with respect of the redshift z G [0,3] 
of the Hubble function, of the matter energy density 
0, and of the deceleration parameter q of the Universe 
filled with magnetic type superconducting dark energy, 
obtained by numerically integrating Eqs. dZni), (El), ©, 
and (1551) . are presented, for a fixed value of cr = 0.0001, 
and for different values of vq, in Figs. [7lfT3l The ini¬ 
tial conditions use to integrate the system of cosmolog¬ 
ical evolution equations are 0(0) = 0o = 0.25, a(0) = 1, 
S(0) = 0.01, 2 ( 0 ) = 5, and S'(0) = 0.01, respectively. 



z 



FIG. 14: Matter energy density of the dust Universe in the 
presence of magnetic type superconducting dark energy as 
a function of the redshift for a = 0.0001, and for different 
values of uq: vo = 0.001 (solid curve), vo = 0.002 (dotted 
curve), Vo ~ 0.003 (short dashed curve), vq = 0.004 (dashed 
curve), and vo = 0.005 (long dashed curve), respectively. 
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FIG. 15: Redshift evolution of the deceleration parameter of 
the dust Universe in the presence of magnetic type supercon¬ 
ducting dark energy for a = 0.0001, and for different values 
of Vq: Vo = 0.001 (solid curve), vo = 0.002 (dotted curve), 
Vo = 0.003 (short dashed curve), vo — 0.004 (dashed curve), 
and Vo = 0.005 (long dashed curve), respectively. 


FIG. 13: The Hubble function of the Universe in the presence 
of magnetic type superconducting dark energy as a function 
of redshift for a — 0.0001, and for different values of vo'- vo = 
0.001 (solid curve), vo — 0.002 (dotted curve), vo = 0.003 
(short dashed curve), vo = 0.004 (dashed curve), and vo = 
0.005 (long dashed curve), respectively. 

As one can see from Fig. [131 the Hubble function of 
the Universe filled with magnetic type superconducting 
dark energy is a monotonically increasing function of 2 
(time decreasing function), indicating an expansionary 
evolution. The matter energy density 0, represented in 
Fig. [T4l monotonically increases with the redshift, and 
tends to zero in the limit of small z. Its dynamics is ba¬ 
sically independent on the adopted numerical values of 
the parameters a and vq. The dust magnetic Universe 


starts from a decelerating state at z = 5, with positive 
values of the deceleration parameter q > 0, shown in 
Fig. |T5| The cosmological evolution is generally decel¬ 
erating for 2 < z < 5, with q reaching the value zero 
at z ~ 2. Then the Universe begins to accelerate, with 
g < 0, and in the large time (small z) limit we have 
limj;_^o'7('2) = —1- Thus, in the final stages of evolution 
of the Universe fillet with magnetic type superconduct¬ 
ing dark energy the cosmological expansion is of de Sitter 
type, with the dark energy driving the Universe’s expan¬ 
sion. The parameter of the dark energy equation of state, 
represented in Fig. 1161 is smaller than zero in the entire 
redshift range 0 < z < 5, and it tends to -1 in the limit 
of small redshifts. 
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FIG. 16: Redshift evolution of the parameter of the dark 
energy equation of state of the dust Universe in the presence 
of magnetic type superconducting dark energy for cr = 0.0001, 
and for different values of uq: uq = 0.001 (solid curve), vo = 
0.002 (dotted curve), vo = 0.003 (short dashed curve), vo = 
0.004 (dashed curve), and vg = 0.005 (long dashed curve), 
respectively. 


IV. THERMODYNAMIC INTERPRETATION 
OF THE SUPERCONDUCTING DARK ENERGY 
MODELS 

In the present Section we analyze the physical inter¬ 
pretation of the superconducting dark energy model by 
adopting the point of view of the thermodynamics of the 
matter creation irreversible processes As we 

have already seen, the energy conservation equation of 
the superconducting dark energy models, Eq. m , con¬ 
tain, as compared to the standard adiabatic conserva¬ 
tion equation, an extra term, which can be interpreted 
thermodynamically as a matter creation rate. According 
to irreversible thermodynamics, matter creation repre¬ 
sents an entropy source, generating an entropy flux, and 
thus modifying the temperature evolution of the consid¬ 
ered gravitational system. On the other hand, due to 
our choice of the geometry of the Universe, all the non¬ 
diagonal components of the total energy-momentum ten¬ 
sor of the superconducting dark energy model are equal 
to zero, so that = Q, ^ In particular, 

from the point of view of the thermodynamics of the ir¬ 
reversible processes and open systems, this condition im¬ 
plies the impossibility of heat transfer in the Friedmann- 
Robertson-Walker models of superconducting dark en¬ 
ergy, since the condition ^ ^ _ 1^2,3 must 

always hold. 


A. Matter creation rates and the creation pressure 


system the second law of thermodynamics, in its most 
general form, is given by 


_d 

dt 


ipa^) 




dQ 

dt 


p + p d 
n dt 


(na^) , 


( 86 ) 


where dQ is the heat received by the system during time 
dt, and n = N/V is the particle number density, respec¬ 
tively. Due to our choice of the geometry of the Universe, 
in a homogeneous and isotropic system filled with super¬ 
conducting dark energy only adiabatic transformations, 
defined by the condition dQ = 0, are possible. There¬ 
fore in the following we ignore proper heat transfer pro¬ 
cesses in the superconductor type cosmological system. 
However, under the assumption of adiabatic transforma¬ 
tions, Eq. (I55|) . representing the general formulation of 
the second law of thermodynamics, contains the term 
[{p+p)/n]d (na^^ /dt, which explicitly takes into account 
the variation of the number of particles in a given vol¬ 
ume. Hence, in the general thermodynamic approach of 
open systems, even for the case of adiabatic transforma¬ 
tions with dQ = 0, there is a ’’heat” (internal energy), 
received/lost by the system, which is entirely due to the 
change in the particle number n. From the cosmological 
perspective of the superconducting dark energy models, 
the change in the particle number is due to the transfer 
of energy from dark energy to matter. Thus in this class 
of cosmological models matter creation acts as a source 
of internal energy, as well as of entropy. For adiabatic 
transformations dQ/dt = 0, Eq. (1861) can be written in 
an equivalent form as 


p + 3(p -I- p)H = ^ ^ (h -I- 2,Hn). (87) 

n 


Therefore, from the point of view of the thermodynam¬ 
ics of open systems, Eq. (EU), giving the energy conserva¬ 
tion equation in the superconducting dark energy models, 
can be interpreted as describing particle creation in an 
homogeneous and isotropic geometry, with the time vari¬ 
ation of the particle number obtained from the equation 


n + 3nH = rn, (88) 

where the particle creation rate T is defined as 


r = (A, - v,<^)] + 

p + p[ 2 

f i [/ (^/5 - V/3<^)] - d^V {A\ </.) - 

2dA2V (A^, (j)) m'^A“V„Aj. \. (89) 


To analyze the thermodynamical implications of the 
superconducting dark energy models at the cosmological 
scale we start with an open system containing N particles 
in a volume V = a^, and characterized by an energy 
density p and a thermodynamic pressure p. For such a 


Therefore the energy conservation equation in the su¬ 
perconducting dark energy model can be written in the 
alternative form 

p + 3{p + p)H = {p + p)T. 


(90) 
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As shown initially in , for adiabatic transformations 
Eq. (1861) . describing irreversible particle creation in an 
open thermodynamic systems, can be formulated as an 
effective energy conservation equation of the form 

4(„„») + (p + pjl„= = 0, (91) 

which can be written in an equivalent form as, 

p + 'i{p + p+Pc)H = (92) 


where we have introduced the terrn^c, called the creation 
pressure, and which is defined as 


Pc 


p + pd[na^) 

n da^ 

p + p r 

3 H' 


P + P 
SnH 


{n + 3nH) = 


(93) 


Therefore in the superconducting dark energy model the 
creation pressure can be obtained as 


f i (^/5 - - 

2dA2V {A^, (f) u''A°‘Vc,A^ \. (94) 


1. Particle creation rates and creation pressure in the 
electric type superconducting dark energy model 

As an example of the thermodynamic description of 
the superconducting dark energy models we consider the 
electric type superconducting dark energy case, for which 
the energy conservation equation is given by Eq. (HD) , can 
be formulated as 


B. Entropy and temperature evolution 


In order to formulate the second law of thermodynam¬ 
ics for open systems, and to apply it to the superconduct¬ 
ing dark energy model, we must decompose the entropy 
change in the cosmological fluid into two components: 
the entropy flow term d^S, and the entropy creation term 
diS. Hence the total entropy S of an open thermody¬ 
namic system can be written as [2^ [2^ 

dS = deS + diS, (98) 

where we assume that diS > 0. Both the entropy flow 
and the entropy production rate in the superconducting 
dark energy model can be evaluated by starting from the 
total differential of the entropy, given by (^ . 

Td {sa^) = d (pa^) + pda^ — pd (na^) , (99) 

where T is the temperature of the open thermodynamic 
system with superconducting particle creation, s = S/a^ 
is the entropy per unit volume, and p is the chemical 
potential, defined in the usual way as 

pn = {p + p) — Ts. (100) 


For closed systems and adiabatic transformations dS = 
0 and diS = 0. However, in the presence of matter cre¬ 
ation there is a non-zero contribution to the entropy. For 
homogeneous systems the entropy flow term deS cancels, 
so that deS = 0. But matter creation from superconduct¬ 
ing dark energy acts as a source for entropy production, 
with the corresponding entropy time variation obtained 

as m 


dt 


^ p + p d 
dt n dt 


(na^) 


S d / r}\ 

T —- (na^) > 0, 
ndt^ ’ - 


(„«=) = 


( 101 ) 


From Eq. (IIOII) we obtain for the time variation of the 
entropy the equation 


a 


p + 3{p + p)H = —p {p -k 3Hp). 


(95) 


dt 


-{ri + 3Hn) = TS>0, 
n 


( 102 ) 


From Eq. (1951) it follows that for p = 0 the matter 
energy is conserved, p + 3{p + p)H = 0, and there is no 
particle creation from the superconducting dark energy. 
However, for p ^ 0, matter and energy transfer processes 
take place in the presence of the superconducting electric 
type dark energy, with the particle creation rate F given 
by 


giving for the entropy increase due to particle creation 
the expression 

= (103) 

where Sq = S{0) is a constant. With the use of Eq. (IMI) . 
we obtain for the entropy production in the supercon¬ 
ducting dark energy models the equation 


T = i-^{p + 3Hp). 


(96) 


4:X p+p 

The creation pressure for this model can be obtained 


as 


=-(III*'*^ 


(97) 


S' dt 


f £ [/ (^/5 - V/3</>)] - d^V (A2, cf) u^A, - 
2dA2V (A^, 4 ,) i > 0. (104) 
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Equivalently, the above equation can be written as 


S~dl 


9 

a p 

4A p + p 


(p + ZHp) > 0. 


(105) 


An important thermodynamic quantity, the entropy 
flux vector of the particles created from the super¬ 
conducting dark energy, is defined according to 


= nau^, 


(106) 


where a = S/N is the specific entropy per particle. The 
entropy flux vector 5^ must satisfy during the entire cos¬ 
mological evolution the second law of thermodynamics, 
which requires that the constraint > 0 be satisfied 

for all times. By taking into account the fundamental 
Gibbs relation |^ . 


thermodynamic variables, the particle number density n, 
and the temperatures T, respectively. Hence the energy 
density p and the thermodynamic pressure p can be ob¬ 
tained, in terms of the particle number n and temper¬ 
ature T, by using the standard form of the equilibrium 
equations of state of the matter created by the supercon¬ 
ducting dark energy, 

p = p{n,T),p = p{n,T). (113) 

Therefore the energy conservation equation Eq. (1^ can 
be written in the expanded form 

^n+^f + 3{p + p)H = {p+p)T. (114) 

With the use of the general thermodynamic relation 


nTda = dp — ^ ^ dn, (107) 

n 

and by using the definition of the chemical potential p of 
the superconducting thermodynamic system as given by 

p = - Ter, (108) 

n 

we obtain 

= {h + 3nH) a -I- = 

^{n + 3Hn) , (109) 

where we have taken into account the important relation 

nTa = p — ^ ^ n = 0, (HO) 

n 

which follows immediately from Eq. (l87)) . With the use 
of Eq. ((M)l we obtain for the entropy production rate due 
to the particle creation processes in the superconducting 
dark energy model the expression 

71 I (y 

= {p + p)t \ + 

f £ [/ - 

2dA^V{A^,4>) 

( 111 ) 

The entropy production rate via superconducting par¬ 
ticle creation processes given by Eq. (Hill) can be written 
in a simpler form as 

In the general case the thermodynamic state of a per¬ 
fect comoving fluid is described by only two essential 


dp p + p T dp 
^ =- (115 

dn n n dl 

we obtain for the temperature evolution of the newly cre¬ 
ated particles in the superconducting dark energy model 
the expression 

l = cl- = cUr-3H). (116) 

1 n 

where the speed of sound Cs is dehned as = dp/dp. 
If the matter newly created from the superconducting 
dark energy satisfies a barotropic equation of state p = 
(7 — 1 ) p, 1 < 7 < 2 , the matter temperature evolution 
can be obtained as 

T = TorC<-^. (117) 

V. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered an electro¬ 
magnetic type dark energy model, in which the elec¬ 
tromagnetic gauge invariance is spontaneously broken. 
The action for such a system must be invariant un¬ 
der gauge transformations, A^(a;) —>■ A^{x) + d^A[x), 
i’nix) —>■ exp {iqnA{x)/h) i/riix). where qn are the charges 
destroyed by the field [i3| • These phase changes lead 
to the formation of an ordered state. By writing all 
charged fields as a function of a scalar field (/{x), when the 
matter fields are integrated out we obtain a Lagrangian 
that is a gauge invariant functional of the fields A^^ and 4 >. 
Such a physical model can explain easily all the observed 
properties of superconductors [la-Ei. Tentatively, we 
also propose it as a dark energy model with a broken 
electromagnetic gauge invariance. 

From a physical point of view the superconducting 
dark energy model is a two-held model, leading to a 
scalar-vector-tensor cosmological theory. It can also be 
viewed as a unihed scalar - vector held dark energy 
model, in which the scalar held (/> and the vector held 
A^ appear in the gauge invariant combination A^ — 
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Moreover, similarly the standard electrodynamic case, we 
have assumed the possible existence of a generalized cou¬ 
pling between the matter current and the gauge invari¬ 
ant combination of the potentials We have 

investigated the cosmological implications of this model, 
by restricting our analysis to the case of a homogeneous 
and isotropic geometry. We have considered two distinct 
classes of models, whose main properties are determined 
by the form of the electromagnetic potential The first 

model corresponds to an electric type choice for the dark 
energy potential, with Aq the only non-zero component. 
For this case the general solution of the gravitational field 
equations was obtained numerically for the dust and the 
radiation filled Universe, respectively. In both cases we 
have assumed that the self-interaction potential of the 
electromagnetic type dark energy is a constant. In both 
cases in the long time limit the Universe ends in a deceler¬ 
ating phase. A similar result is obtained for the magnetic 
type dark energy model, in which the electromagnetic po¬ 
tential is restricted to the form = (0, A(t), A(t), A{t)). 
For this case we have investigated, by numerically solving 
the gravitational field equations, the zero thermodynamic 
pressure (dust) cosmological model, in the presence of a 
constant self-interaction potential. Similarly to the elec¬ 
tric type dark energy model, the magnetic superconduct¬ 
ing dark energy model drives the Universe, in the long 
time limit, in an accelerated, de Sitter type, expansionary 
phase. 

Due to the coupling between the dark energy potentials 
and the matter current, the matter energy-momentum 
tensor is not conserve in the present approach. We have 
interpreted, in the framework of the thermodynamics of 
open systems and irreversible processes the non¬ 

conservation of the matter energy-momentum tensor as 
describing particle creation, and energy transfer from the 
superconducting dark energy to ordinary matter. We 
have explicitly obtained the particle creation rates, as 
well as the effective creation pressure generated by the 
irreversible transformation of the field energy into mat¬ 
ter. The entropy production rate and the overall entropy 
evolution of the Universe was also obtained, with the to¬ 
tal entropy being given by the exponential of the time 
integral of the particle production rate F. 

The possible observational study in a cosmological con¬ 
text of the irreversible matter-creation processes in the 
homogeneous and isotropic flat Friedmann-Robertson- 
Walker geometry in the superconducting dark energy 
models may represent one of the possibilities of consider¬ 


ing the viability of this dark energy model. However, in 
order to confirm the validity of the superconducting dark 
energy model developed in the present paper, it is neces¬ 
sary to carefully consider a much wider range of cosmo¬ 
logical and astrophysical tests for this type of models. In 
particular, an essential test of the superconducting dark 
energy model would be the investigation of the classi¬ 
cal macroscopic predictions of the model in large-scale 
structure formation with linear perturbations, and with 
the consideration of the Newtonian limit for small scales. 
Supernovae observations fitting, and the study of the ef¬ 
fects of the matter creation on the Cosmic Microwave 
Background anisotropies could lead to other important 
tests and parameter constraints of the model. Essen¬ 
tially the superconducting dark energy model introduced 
in the present paper is a simple toy model, whose main 
goal is to stimulate the study of alternative or more gen¬ 
eral electromagnetic type dark energy models. 

The superconducting dark energy model introduced in 
the present paper leads to the possibility that matter cre¬ 
ation, associated with matter current - dark energy elec¬ 
tromagnetic potential coupling may also happen in the 
present - day universe, as initially considered by Dirac 
[28l |. The existence of some forms of coupling between 
matter and dark energy are not in contradiction with 
the cosmological observations or with some astrophysical 
data [ 2 ^. However, firm observational evidence for parti¬ 
cle creation on a cosmological scale is still missing. Hope¬ 
fully, a key ingredient of the present model, the functional 
form of the self-interaction potential V (, (()), which 
essentially determines the cosmological dynamics in the 
superconducting dark energy model, will be provided by 
fundamental particle physics (or perhaps even condensed 
matter) models, thus permitting an in depth comparison 
of the predictions of the model with high precision ob¬ 
servational cosmological and astrophysical data. 
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Appendix A: The divergence of the 
energy-momentum tensor 

In order to obtain the divergence of the matter energy- 
momentum tensor in the superconducting dark energy 
model we compute first the divergence of the electromag¬ 
netic type term, 

^ ^ . (Al) 

Hence for the divergence of the electromagnetic com¬ 
ponent we find first 

(A2) ^ 

By taking into account that (I/Jtt) = J“, and 

^vFap = — Vaf/3i/ — y it follows that 

— ^^1 2 pv 2 ^ p^va 
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pP-a\7 p \ — F T° 

^ V UOt I ^ vet •J 


(A3) 


The terms in the bracket vanish, and therefore we find 


V _ p 

V y - ^ 


{Ap - X/p^) + - 


a 


XF^, (A^ - V^cj^) , 


(A6) 


and 


otj^ (V^Aj, - XJfjXJ^cjX) - {XJ^Ap - XJ^XJpej)) = 
(V^A, - V^V.<^) + (A7) 


2dA2V(A^,<f)) A“ 


(A4) 


Then for the divergence of the matter energy- 
momentum tensor we obtain 


- F,^r F X (V^A'^ - X/^X/^cj,) (A, - X/^cj,) + 

A (A'^ - V^(/.) F^, + aS/^f (A, - V,(/.) + 

(V^A. - X/^X/.4>) + ^fF^, - |V./ [Ap - X/p4>) + 

dpV (A2, cf) X7,cp + 2A<^dpV (A2, <).) V,A„ = 0, (A5) 


respectively. 

With the use of the evolution equation of the scalar 
field we find the relation, 

A (V^A'^ - V^V^<(.) (A, - V.(/)) = 

d^V (A2, </,) (A, - (A, - . (A8) 

By substituting the above relation and the expression of 
the divergence of the electromagnetic part of the energy- 
momentum tensor in Eq. (|A8I) . we finally obtain 


where we have used the identities 

A (A^ - {X7^A, - X7^X7,^) - ^ (V,A“ - V,V“0) x 

(A„ - V„<^) - ^ (Al“ - V“0) (V,A« - = 


ir (A, - V,<^)] - {Ap - Vpej,) + 

d^V {A^,(t>) A, + 2dA2V (A2,0) A“V„A, = 0. (A9) 




